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A microscopic study of the individual annealed (In,Ga)As/GaAs quantum dots is done by means of high- 
resolution transmission electron microscopy. The Cauchy-Green strain-tensor component distribution and the 
chemical composition of the (In,Ga)As alloy are extracted from the microscopy images. The image process¬ 
ing allows for the reconstruction of the strain-induced electric-field gradients at the individual atomic columns 
extracting thereby the magnitude and asymmetry parameter of the nuclear quadrupole interaction. Nuclear mag¬ 
netic resonance absorption spectra are analyzed for parallel and transverse mutual orientations of the electric- 
field gradient and a static magnetic field. 

PACS numbers: 78.67.Hc, 78.47.jd, 76.70.Hb, 73.21.La 


INTRODUCTION 

The spin physics of semiconductors has been developed 
for bulk materials and demonstrated a wide variety of lin¬ 
ear and nonlinear phenomena, realized thanks to the optical 
orientation.^ It has been reborn in semiconductor nanostruc¬ 
tures in recent decades.^ A significant part of the spin-related 
phenomena are underlain on the dynamic spin polarization of 
the nuclear spins being polarized by means of the transfer of 
the photon angular momentum to the nuclear-spin system via 
electron-nuclear hyperfine interaction.^’"^ The achievement of 
a high polarization of the nuclear-spin system becomes chal¬ 
lenging for quantum dot (QD) systems where a single spin of 
an electron would be strongly localized^ and is under the in¬ 
fluence of the nuclear spin fluctuations paving the way to a 
fast carrier-spin relaxation.^'^^ In spite of a combination of a 
large variety of methods tried to be used for reaching a suffi¬ 
ciently high degree of nuclear-spin order, the experimen¬ 

tal achievement of the spin polarization, close to hundred per¬ 
cent, is still a challenging problem, limited, in some particular 
cases, by the quantum nature of the spin system.A fur¬ 
ther microscopic analysis including a combination of experi¬ 
mental methods like nuclear magnetic resonance (NMR)^^“^^ 
and/or various spin noise measurements^^""^^ in combination 
with atomistic modeling"^ are highly required to gain in¬ 
sight into such an intriguing problem. 

In self-assembled QDs, A ~ 10^ nuclear spins interact with 
the localized-electron spin with different strengths because of 
a spread of the electron density. However, an additional spread 
of the interaction exists even at nuclear-spin level due to the 
crystal lattice deformation caused by built-in strain. This in¬ 
troduces a non-homogeneous nuclear quadrupole interaction. 


changing the usual nuclear-spin dynamics"^^ but also modify¬ 
ing the NMR spectrum both in the single QDs^^ and in the QD 
ensembles.The microscopic analysis of the quadrupole 
interaction is a rather complex problem while any direct ex¬ 
perimental measurement of its magnitude would hardly be re¬ 
alized in practice. Since only scant experimental progress can 
be expected thereupon, an atomistic analysis would at least 
give some microscopic information within the framework of a 
chosen model.'^^ 

In this paper, we address the investigation of the struc¬ 
tural properties of a single QD with respect to the analysis 
of nuclear quadrupole interaction. We make use of the high- 
angle annular dark-field (HAADF) imaging of an individual 
self-assembled QD in a scanning transmission electron micro¬ 
scope (STEM). This allows us to resolve a crystal lattice with 
atomic-column resolution that can be used to determine the 
shape and chemical composition of the QD. The investigated 
sample is a heterostructure grown by molecular-beam epitaxy 
on a GaAs substrate and contains 20 layers of (In,Ga)As QDs, 
embedded in a GaAs matrix. The post-growth thermal an¬ 
nealing of the structure allows for activation of in-diffusion 
of Ga atoms inside the QD that reduces the number of struc¬ 
ture defects and, in addition, blue shifts the ground-state ex- 
citonic transition.^^""^^ The obtained STEM images were post- 
processed with geometric phase analysis (GPA), from which 
the Cauchy-Green strain tensor components are extracted. In 
addition, to obtain the profile of the In and Ga concentra¬ 
tions inside the QD, the local chemical composition is de¬ 
termined by means of atomically resolved HAADE-STEM 
and energy-dispersive X-ray (EDX) spectroscopy. These data 
were then used to reconstruct the distribution of the strain- 
induced electric-field gradients (EEG) causing the nuclear 
quadrupole interaction. Eurther calculations allow the inter- 
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Figure 1. Cross-section micrograph of the heterostructure with an¬ 
nealed (In,Ga)As/GaAs QDs obtained using scanning transmission 
electron microscopy. 


pretation of the NMR transitions that can be detected in the 
ensemble of such 

The paper is organized as follows. In Sec. I, we present the 
details of the microscopic study. The results of the postpro¬ 
cessing of the microscopic images are discussed. In Sec. II, 
we discuss the possibility of application of microscopic meth¬ 
ods for an analysis of strain-induced nuclear quadrupole in¬ 
teraction and quantitative modelling of the NMR spectra. The 
main results are summarized in Sec. III. 


I. HAADF-STEM IMAGING AND STRAIN MAPPING 

To obtain the crystal strain and, correspondingly, the 
quadrupole interaction, high-quality microscopic imaging is 
required.^^ In particular, the image should include the inves¬ 
tigated QD and a region of the surrounding GaAs matrix of 
much larger area than the size of the QD. This region is used as 
a bulk reference area for further analysis. Selection of the ap¬ 
propriate microscopy method is based on the ability to resolve 
the atomic columns of the crystal and the ability to extract 
the chemical composition of the QD. Such data are obtained 
by analysing the heterostructure with STEM using HAADF- 
detector, providing information about the concentration of the 
atoms of different species. 

The cross-sectional specimens of the sample with an area of 
about 0.2 cm^ are prepared by mechanical cutting of the het¬ 
erostructure in the (lIO) crystallographic plane and then glue¬ 
ing to a tripod holder. The specimens are then mechanically 
polished reducing the thickness of the sample down to ~10- 
30 [im. Further ion etching of the QD containing region is ap¬ 
plied using a Gatan precision ion polishing system (PIPS) un¬ 
til a hole is etched into the sample. At the edge of this hole the 
sample thickness is clearly below 200 nm and thus transpar¬ 
ent for electrons. The experiments are performed with an FEI 


Figure 2. High-resolution STEM image of a single (In,Ga)As QD 
embedded in the GaAs matrix. The shape of the QD is marked by the 
thin solid line guide to the eyes with respect to the Indium concen¬ 
tration profile shown in Fig. 3. Atomic-column resolution is demon¬ 
strated in the inset. 


Titan 80/300 microscope operated at 300 kV electron-beam 
acceleration voltage and equipped with an image-aberration 
corrector, EDX, and HAADF detectors. The overview image 
of the structure is shown in Fig. 1, from which multiple QDs 
placed on several wetting layers are resolvable. Size and shape 
of the dots are found to be weakly dispersed, which is also 
confirmed by a relatively narrow linewidth of the ground-state 
photoluminescence studied previously with these QDs."^^’^^ A 
few randomly chosen single QDs are analyzed. 

Consequently, the high-resolution HAADF-STEM images 
were recorded in the QD-containing region of the sample. The 
image size is selected to approximately 50 nm x 50 nm and the 
spatial resolution of the microscope is around 0.12 nm in the 
scanning mode. Each image contains one QD plus a sufficient 
amount of the surrounding GaAs matrix, which can be used 
as reference material for thickness measurement. As high 
resolution was obtained in the HAADF-STEM images, they 
were not only suitable for atomic Z-contrast evaluation^^"^^ 
but could be also used for geometric phase analysis.^^’^^ The 
chemically sensitive image contrast of the HAADF-STEM 
images (see Fig. 2) is due to the used HAADF-detector. This 
ring-shaped detector detects only electrons that are scattered 
into high-angles (36-220 mrads for the camera length used 
in this work). The amount of electrons scattered into this re¬ 
gion strongly depends on the nuclear charges of the scattering 
specimen atoms. 

The following evaluation procedure was performed for ev¬ 
ery single image in order to obtain the chemical composition. 
First, all atomic columns in the high-resolution image were 
identified. To this end, a Wiener Filter has been applied on 
the original image."^^’^^ After that, the image has been divided 
into “Voronoi”-cells and the mean intensity has been calcu¬ 
lated for each “Voronoi”-cell and has been assigned to the 
corresponding atomic column. Note that in this step the inten- 
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Figure 3. (Color online) Indium concentration map in (In,Ga)As 
solid compound evaluated from HAADF-STEM. The solid line rep¬ 
resents a boundary of the QD surrounding the piece of sample with 
an In-concentration higher than 0.1 reaching the maximum value 0.3 
in the center of the QD. 


sides of the original unfiltered image were used. To allow for 
the quantitative comparison to simulated data, the mean inten¬ 
sity values were normalized with respect to the intensity of the 
scanning electron probe. In the following step, the normalized 
intensities were compared to reference data from multislice 
simulations in the frozen lattice approach using the STEMsim 
program^^ carried out in dependence of sample thickness and 
indium fraction. For more details on the simulations we re¬ 
fer the reader to Ref. 57. Finally, the sample thickness was 
evaluated from the GaAs region and interpolated over the QD 
to allow for the determination of the indium concentration as 
shown in Fig. 3. 

To verify the accuracy of HAADF-STEM, EDX spec¬ 
troscopy is used. The EDX spectra were also acquired in 
HAADF-STEM mode. While the beam is centred on a small 
region of the sample surface, the X-ray counts were integrated 
over one minute. Both methods of the Z-contrast evaluation 
show very similar values for the maximum concentration of 
Indium in the single QD, i. e., 0.3 + 0.05 and 0.35 + 0.05 for 
EDX and HAADF-STEM, respectively. 

Several methods for mapping the crystal-lattice strain can 
be applied.^^’^^’^^"^^ In the electron microscope, none of 
them is able to evaluate the strain-tensor component along 
the electron-beam propagation direction. However, the two- 
dimensional maps of two diagonal strain-tensor components 
could be extracted. For further analysis, the component of 
shear strain is also required. Denoting the Cartesian coor¬ 
dinate system with respect to crystallographic axes so that 
X [110] and z [001], the Cauchy-Green strain-tensor 
components and s^z are obtained. While other com¬ 

ponents are experimentally unaccessible, the symmetry of the 
problem allows us to equalize the x and y directions, keep¬ 
ing in mind that the specimen is several times thicker than 


the QD lateral size. We take Syy = s^x and Sy^ = Sxz keep¬ 
ing £xy - 0. While Exy = 0 is a reasonable assumption since 
it is large at the heteroboundary, which is unsharp in the an¬ 
nealed QDs, other assumptions emerge as a value judgement. 
Note, that in (In,Ga)As/GaAs QDs, there is a small inequiv¬ 
alence of [110] and [lIO] directions revealing in the exciton 
fine-structure splitting"^^ that can be changed in strained QDs, 
as advanced atomistic modeling shows. 

The strain mapping is performed by means of two methods, 
the geometric-phase analysis (GPA)^^ and a modified peak- 
pair analysis (PPA)^^ for verification. The CPA allows us 
to extract the crystal-lattice distortions locally at each atomic 
column with respect to the unstrained lattice area. The basic 
idea of such an analysis is based on the fast Fourier transform 
of the real-space image into the reciprocal space. Focal strain 
tensor components, symmetric and rotation parts of distortion 
are calculated by derivation of the displacement obtained from 
two non-collinear Fourier components.^"^ The phase compo¬ 
nent of this function, called the geometric phase, describes 
the position-dependent lattice deviation with respect to a ref¬ 
erence. The reference area of the crystal lattice is taken in the 
barrier region 40 nm away of the QD where the In concentra¬ 
tion is negligibly small. The CPA strain-tensor components 
are extracted as 


= -^G-iV(D,(r), (1) 

where G = [ matrix formed by the components 

of two non-collinear reflexes gi and g 2 , each of which is con¬ 
nected with the position-dependent geometric phase, 0 ^(r), 
as 0^(r) = 2;rAg(r) • r, where g(r) represents the periodicities 
corresponding to the Bragg refiections. It has the following 
relationship with the lattices fringe spacing: J = l/|g|. 

Additionally, the same STEM image with a single QD is 
processed using the PPA. In contrast to the CPA, the PPA is 
a real space procedure for strain mapping. PPA works with 
images having well-resolved fringe patterns, finding pairs 
of peaks along a preselected direction and distance in the 
affine transformed space defined by a pair of basis vectors 
a = (Gx, Gy) and b = {bx, by). When the reference vectors 
are chosen on the filtered image, they can be used to construct 
an affine transformation. The next step in the PPA is the iden¬ 
tification of peak-pairs using the chosen basis vectors and the 
intensity maxima set in the image. The Cauchy-Green strain 
components can be calculated by solving the following set of 
linear equations 
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with coordinates of the displacements {Ux, Uy) and (Vj^, V 3 ;) with 
respect to the reference vectors a and b. 

The methods described above are based on the same gen¬ 
eral assumption that the relation between the atomic column 
positions in real crystal and features in the STEM image is 
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Figure 4. (Color online) Cauchy-Green strain tensor components extracted from the GPA analysis of HAADF-STEM images of a single QD, 
from left to right: (a) (b) and (c) s^z- 


constant inside the studied region, where the phase shift be¬ 
tween maxima and atoms is supposed to be constant. Both of 
these two methods, within the calculation error, gave identi¬ 
cal results for all three tensor components Sxx, and s^z for 
the single QDs under study. The PPA requires less amount of 
memory and calculation time, given that the two-dimensional 
complex Fourier transform is not required. On the other hand, 
PPA fails when lattice peaks are not easily detected due to 
lower resolution of the image and appearance of sublattices 
due to a structure defect. In this case, filtering of the origi¬ 
nal image can be performed. Summarizing, we consider that 
both algorithms, GPA and PPA are useful for strain mapping, 
each having different advantages and limitations, and should 
be considered in each particular case specifically. In this work, 
we present the results of GPA only. 

In general, the strain is measured with respect to the bulk 
lattice parameters of the material. Following this definition, 
the GPA (PPA) strain is connected with the material strain, Sij, 
as follows 


, GPA(PPA)^ A (^ / 

(l+e,. V)) = (l-%(r))—(3) 


Here, the position-dependent lattice constant is determined by 
taking into account the concentration dependence of the solid 
compound via Vegard’s law, ^(in,Ga)As^j.^ _ ^GaAs. _ ^(r)) + 

^inAs. where c(r) is the position-dependent In concentra¬ 
tion shown in Fig. 3. 

Figure 4 illustrates the components of the physical strain, 
and Sxz evaluated from GPA analysis. The strain- 
tensor components extracted from high-resolution STEM im¬ 
ages are averaged over the sample thickness, i.e., their values 
in the QD would be slightly larger in magnitude. Note, that an 
additional source of strain release arises from the sample thin¬ 
ning down to less than a qm. In this case, the QD becomes 
closer to a surface than a deep-in-bulk QD. Therefore, both 
these effects should be considered as a source of measurement 
error for further evaluation of the results. 

The negative value of Sxx shows that the crystal lattice is 
compressed in x direction ([110]) [see fig. 4 (a)]. On the 
contrary, the s^z. component is mostly positive in the upper 


left region of the QD (see fig. 4 (b)), i.e., there is a stretch¬ 
ing of the lattice in the growth direction z ([001]). Fig. 4(c) 
shows a shear-strain Sxz, which is present in a relatively small 
area having a mostly positive value and some periodicity in 
its distribution. In the wetting layer, the values of shear strain 
are predominantly negative, while inside the QDs this defor¬ 
mation is positive. Overall, the obtained results qualitatively 
well coincide with the finite-element modelling, represented 
previously in Ref. 45. 


II. ANALYSIS OF STRAIN-INDUCED QUADRUPOLE 
INTERACTION 

The change of energy experienced by the nuclear spin I can 
be dependent of nuclear orientation. The charge environment 
of the nuclei interacts with the external electric potential V. 



I—I 0.0 
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Figure 5. (Color online) EFG averaged over the propagation direc¬ 
tion of the electron beam. The red arrows indicate the relative value 
of the strain-induced quadrupole frequency. All the microscopy data 
errors are indicated by the additional black arrows. The surface plot 
shows the spatial distribution of the EEG asymmetry 77 . Black lines 
indicate the shape of the QD extracted from the (In,Ga)As alloy con¬ 
centration map. 
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Figure 6 . (Color online) The field-frequency maps of the nuclear spin-flip transition rates given by Eq. ( 6 ) calculated when the magnetic field 
is oriented perpendicular (a) or parallel (b) to the QD growth axis. The raster map indicates the probability of the transitions for ^^Ga nuclei 
with spin 7 = 3/2 (as an example). The transitions between dipole-allowed and dipole forbidden states calculated with 77 and Vzz the QD 
volume are marked by solid and dashed lines, respectively. The inset graphs demonstrate the splittings of the nuclear-spin states for a single 
nucleus affected by a mean vq in symmetric (dashed lines) or asymmetric (solid lines) quadrupole configurations. 


In the equilibrium, the nuclei experience zero average electric 
field but nuclei having 7 > 1 /2 have also nonzero quadrupole 
moment interacting with the EFG denoted hereafter as 
The physical origin of nonzero Vtj is any inhomogeneity of the 
electric fields. The strongest effect is caused by the substitu¬ 
tional atoms in the alloy. Due to different electronegativity of 
the cation atoms, the anion As nuclei are highly affected by 
the EFG caused by the charge environment redistribution. For 
example, 5% electronegativity difference between the group- 
ill atoms leads to the nuclear-spin splitting of several tens of 
MHz, as shown in (Al,Ga)As bulk.^ 

In self-assembled QDs, a weaker splitting would be ob¬ 
tained due to the crystal-lattice biaxial strain of several per¬ 
cent, leading to the quadrupole splitting of nuclear-spin states 
of the order of a MHz."^^ Even in the absence of both effects, a 
strong electron localization would lead to the additional inho¬ 
mogeneity of spin splitting caused by interaction with the ap¬ 
plied electric field Our estimation shows, however, that 
the EFG variation due to inhomogeneity of the ground-state 
electron density inside the model QD leads to, at least, an or¬ 
der of magnitude smaller variation of spin splitting than its 
strain-induced value. However, in electrically driven QDs or 
in self-assembled QDs with a bias-controlled charge state, the 
electric-field induced EFG variation requires additional veri¬ 
fication, as shown previously in a quantum well system.^ ^ 
Considering further the strain-induced quadrupole interac¬ 
tion, the Hamiltonian reads as 

= -hyA ■ B + ^( 3 /^ - r + + /^)). ( 4 ) 

Here, the first term couples nuclear spin I of the 7th nucleus 
having the gyromagnetic ratio hyj with the external magnetic 
field B. The second term describes the nuclear quadrupole in¬ 
teraction with strength hvQ. The coordinate system (X, J/, Z) 
corresponds to the major frame of the EFG acting on the nu¬ 
cleus generally not coinciding with the (x, y, z) frame of the 
crystal. The quadrupole frequency vq is proportional to the 


major EFG component Vzz satisfying Vzz > 

The quadrupole asymmetry parameter, 77 = (Vxx-yyy)/yzz^ 
is determined by the relative ratio of its other principal com¬ 
ponents. 

The microscopic analysis allows to quantitatively evaluate 
the EFG tensor and, as a consequence, to map the quadrupole 
interaction over the sketch of the QD. The EFG tensor re¬ 
lates to the elastic strain as Vtj = Sijkm^km^ where Sijkm 
is the fourth-rank gradient elastic tensor.Using the strain- 
tensor components extracted from HAADF-STEM analysis 
(see Fig. 4), we calculate all non-zero components of Vtj, and, 
correspondingly 


hvQ = 


eQVzz 

4/(2/- 1) 


(5) 


where Q is the nuclear quadrupole moment. The results of 
such an analysis for ^^Ga are plotted in Fig. 5. The magnitude 
of Vq varies by the order of magnitude over the QD volume. 
However, the principal direction Z. of the EFG varies by ap¬ 
proximately 10 degrees, keeping a mean value misalignment 
of 8 degrees from the [001] crystal axis (see Fig. 5). The vari¬ 
ation of the EFG asymmetry 77 is found to be less than 0.2 over 
the whole QD volume. 

To quantitatively evaluate the NMR absorption, the eigen¬ 
value decomposition of Eq. (4), taking into account the spatial 
variation of vq and 77 , is made. The spin-flip rates between |/) 
and 1 7 ) eigenstates having energies Ei and Ej are expressed as 
follows^^ 


2jTT O 

= z ( 6 ) 


when the z axis coincides with the quantization axis of 
electron-spin observable. The transitions can be induced 
by either the absorption of radio-frequency electromagnetic 
fields or by the temporal variation of the Knight field due 
to electron-spin pumping, both modulated at frequency /mod- 
The eigenstates |/) and I 7 ) are no more eigenstates of /, there¬ 
fore the transitions with momentum projection by Attz/ > 1 
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are allowed under certain conditions'^ The transition rates 
are calculated when scanning the magnitude of the external 
magnetic field oriented perpendicular [Fig. 6 (a)] or parallel 
[Fig. 6 (b)] to the [001] crystallographic axis, respectively. The 
statistics caused by the inhomogeneous variation of vq and 77 
shown in Fig. 5 over the QD is taken into account here. 

In both orientations of the magnetic field, the transition 
lines are strongly broadened due to inhomogeneity of the 
quadrupole interaction within the QD volume. A geometry 
with B aligned across the growth axis is more sensitive to the 
value of vq. Here, the Kramers doublets with m/ = +1/2, 
±3/2, etc., instantaneously mix in the B field, thus splitting 
the states linearly for m/ = +1/2 and nonlinearly for higher 
states. The additional asymmetry of the quadrupole config¬ 
uration given with 77 > 0 or a small tilt of the magnetic field 
from exactly perpendicular to the EFG axis results in a change 
of the splitting, particularly for states having \mi\ > 1/2. 
Both, the variation of vq and tj, and the magnetic-field tilt¬ 
ing result in a broadening of the transition lines, as Fig. 6(a) 
shows. As a characteristic scale, the strongest transition be¬ 
tween 777/ = +l/2of^^Ga spreads over about several hundred 
kHz in moderate fields B > 40 mT [see Fig. 6(a)]. If the 
magnetic field is oriented along the EFG axis, the nuclear- 
spin states are eigenstates of the Hamiltonian and do not mix 
by the magnetic field. However, a small tilt of the field re¬ 
sults in a mixing, causing an anticrossing of lowest-energy Ei 
and E 2 states [see inset in Fig. 6(b)]. Further mixing of states 
is provided by the asymmetric quadrupole configuration with 
77 > 0. Both effects lead to broadening of the sharp lines to 
values, however, smaller than in the perpendicular geometry 
by a factor of two, at least. This goes well with the recent 
observation of NMR lines in a single QD reported in Ref. 42 
where the tail of the line of several hundred kHz is observed 
for all nuclear isotopes. 

The shear strain expressed in the dot as providing 

nonzero rj further modifies qualitatively the NMR spectra. For 
the sake of comparison, the energy-level schemes in both sym¬ 
metric and asymmetric quadrupole configurations are shown 
in the insets of Figs. 6 (a) and 6 (b). The frequency shifts of 
several tens of kHz are observed in both geometries. In addi¬ 
tion to that, the spin-level E 2 and £"3 anticrossing exists when 
77 > 0 [see Fig. 6 (b)]. The sample tilting is a natural way to re¬ 
duce the broadening,^^ that works well in single-QD NMR."^^ 
However, full compensation of the spin-level anticrossings is 


prohibited if 77 9 ^ 0, as demonstrated in Fig. 6 (b). 

III. CONCLUSION 

To conclude, the microscopic study of the quadrupole inter¬ 
action in self-assembled (In,Ga)As QD is done by using scan¬ 
ning transmission electron microscopy. The HAADF-STEM 
technique allows us to extract the (In,Ga)As alloy concentra¬ 
tion of In and Ga atoms, keeping the atomic-column resolu¬ 
tion of the microscopy image, and to evaluate the in-plane 
components of the Cauchy-Green strain tensor by using the 
geometric phase analysis. Further mapping of the biaxial and 
shear strain components allows for quantitative reconstruction 
of the strain-induced EFG tensor components. Modelling the 
NMR absorption spectra, the magnitudes of the NMR lines 
broadenings and a shift of certain NMR transitions caused by 
the asymmetry of the EFG tensor are evaluated. In particular, 
the asymmetry parameter of the quadrupole interaction leads 
to shifts and anticrossings of certain transitions in the NMR 
absorption spectrum. 
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